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Abstract 

The Periodic Anderson Model (PAM) in the U — ► 00 limit has been studied in a 
previous work employing the cumulant expansion with the hybridization as pertur- 
bation (M. S. Figueira, M. E. Foglio and G. G. Martinez, Phys. Rev. B 50, 17933 
(1994)). When the total number of electrons Nt is calculated as a function of the 
chemical potential /x in the "Chain Approximation" (CHA), there are three values 
of the chemical potential [i for each Nt in a small interval of Nt at low T (M. S 
Figueira, M. E Foglio, Physica A 208 (1994)). We have recently introduced the "X- 
boson" method, inspired in the slave boson technique of Coleman, that solves the 
problem of non conservation of probability (completeness) in the CHA as well as 
removing the spurious phase transitions that appear with the slave boson method 
in the mean field approximation. In the present paper we show that the X-boson 
method solves also the problem of the multiple roots of Nt(/x) that appear in the 
CHA. 
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1 Introduction 



In the Periodic Anderson Model (PAM) there are four local states at each site 
j of the lattice (identified by a single index in this work): the vacuum state 
| j, 0), the two states | j, a) of one electron with spin component a = ±1 and 
the state | j, 2) with two local electrons [1]. In the limit of infinite Coulomb 
repulsion U — > oo, the state | j, 2) is empty, and we shall use the Hubbard 
operators [2] to project it out from the space of local states at site j. The 
Hamiltonian of the system is then [3]: 



k,cr j,a 
E O^.k-VU'-k,, + V^claXjfia) , (1) 

i,cr,k 

where the first term is the Hamiltonian of the conduction electrons (c-electrons) , 
the second term describes independent localized electrons (/-electrons), and 
the last term is the hybridization Hamiltonian giving the interaction between 
the c-electrons and the /-electrons. 

The X Hubbard operators do not satisfy the usual commutation relations so 
that diagrammatic methods based on Wick's theorem are not applicable, and 
one has to use product rules instead: 

Xj ia b-Xj iC d = 5b,cXj : ad- (2) 

The identity decomposition in the reduced space of local states at site j is 
then 

Aj,oo + X jta<7 + Xj^ = Ij, (3) 

where a = —cr, and the three X^ aa are the projectors into | j,a). Because of 
the translational invariance, the occupation numbers n^ a =< X^ aa > satisfy 
n,j ta =n a (independent of j), and from Eq. (3) we obtain the "completeness" 
relation 

n + n a + n w = 1. (4) 

The occupation numbers can be calculated from appropriate Green's func- 
tions (GF), and it has been found that Eq. (4) is not usually satisfied when 
the n a are calculated with approximate cumulant Green's functions (GF) [4]. 
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An approximation with this behavior is the "Chain Approximation" (CHA), 
which was first obtained by Hewson [5,6]. This approximation has interesting 
properties: it is $— derivable [7,8], and it is also the most general cumulant 
expansion with only second order cumulants. We have employed several proce- 
dures to restore completeness to the CHA: renormalization of the one-electron 
Green's functions (GF) or adding diagrams to the GF in the CHA. The second 
method lead us to a conjecture on a systematic way of achieving complete- 
ness by adding a set of diagrams to an arbitrary family [8] . An alternative to 
these techniques was inspired in the mean field treatment of the slave boson 
technique [9,10], in which the correlated problem is transformed into an un- 
corrected one with one condition that forces to zero the occupation of | j, 2) 
states. This condition turns out to be just our Eq. (4), and following Coleman 
we minimized the free energy of the system calculated with the CHA but forc- 
ing the validity of Eq. (3), that implies the completeness. This is the essence of 
the X-boson method [11,12] and, differently from the slave boson treatment, 
the correlation is kept at the final stage because it is intrinsic to the CHA [13]. 
An important consequence of this fact is that the spurious phase transition 
that appears in the slave boson treatment for several regions of the system 
parameters (for intermediate temperatures or when fi» Ef) [14], disappears 
completely from our treatment. The results of the X-boson method are fairly 
close to those obtained by the slave boson treatment in the region of its va- 
lidity, while the X-boson method gives results that are close to those of the 
CHA when the slave boson is not valid any more (at high temperatures and 
when fx » ej ). 

A rather inconvenient aspect of the CHA is that its results show a region of 
instability [15], apparent because the dependence of N t with fi shows multiple 
values of fx for a given N t within a small interval of values of Nt, and this 
leads to negative compressibility. The main result of the present work is to 
show that this difficulty of the CHA is removed by the X-boson treatment, 
and N t (fx) becomes then a monotonous function. 

The paper is organized as follows: the X-boson approach is presented in section 
2, the GF in the CHA are given in section 3, and in section 4 we describe the 
calculation of the free energy. Our results and conclusions are presented in 
section 5. 



2 X-boson Cumulant Method 

In Coleman's "slave boson" method [9,10], the Hubbard X operators are writ- 
ten as a product of ordinary bosons and fermions: Xj,oo bjbj, Xj )0a — > 
bj~fj tC7 , Xj tUO — > f^bj, and a condition, that is equivalent to Eq. (4), is im- 
posed to avoid states with two electrons at each site j. In the spirit of the 
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mean field approximation bf ^< bf >= r and the method of Lagrangian 
multipliers is employed to minimize the free energy subject to that condition. 
The problem is then reduced to an uncorrelated Anderson lattice with renor- 
malized hybridization V — > rV and / level £/—>£/ + A, and the conservation 
of probability in the space of local states is automatically satisfied because 
they are described by Fermi operators. 

The approximate GF obtained by the cumulant expansion [3,16] do not usually 
conserve probability (i.e. they do not satisfy Eq. (4)), and the procedure we 
adopt to recover this property in the X-boson method is to introduce 

R =< X jyOQ >=< bpj >, (5) 

as variational parameter, and to modify the approximate GF so that it mini- 
mizes an adequate thermodynamic potential while being forced to satisfy Eq. 
(4). To this purpose we add to Eq. (1) the product of each Eq. (4) into a La- 
grange multiplier Aj, and employ this new Hamiltonian to generate the func- 
tional that shall be minimized by employing Lagrange's method. To simplify 
the calculations we use a constant hybridization V, as well as site independent 
local energies £/j> = and Lagrange parameters Aj = A. We then obtain 
a new Hamiltonian with the same form of Eq. (1): 

H = Y E t ct cr + 

J2ef,*Xj,v<r + N s A(R-l) + 

v E (4°. + C L , (6) 

but with renormalized localized energies 

e f ,a = £f,c + A. (7) 

The parameter 

R = 1 - £ < X CTCT > (8) 

CT 

is now varied independently to minimize the thermodynamic potential, choos- 
ing A so that Eq. (4) be satisfied. While at this stage the electrons in the 
slave boson Hamiltonian have lost all the correlations, the Eq. (6) is still in 
the projected space and it is not necessary to force the correlations with an 
extra condition. On the other hand we do not have an exact solution for this 
new problem, and we then consider the most simple approximation obtained 
within the cumulant formalism, the Chain approximation (CHA) [5,6]. The 
need of minimizing a thermodynamic potential arises because the complete- 
ness relation is not automatically satisfied for approximate cumulant solutions, 
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and although the two procedures are formally very similar, they have a rather 
different meaning. 

The Grand Canonical Ensemble of electrons is employed in the present treat- 
ment, and instead of Eq. (f ) we have to use 

H-H-AY, cl,aC k , a + £ u a X haa \ , (9) 

I k,cr ja ) 

where v a — 0, 1 is the number of electrons in state | a >. It is then convenient 
to define 

£ j,a = E f,j,a - (10) 

and 

^kcr = E ko - — H , (f 1) 

because -E/,j, a and appear only in that form in all the calculations. The 
exact and unperturbed averages of the operator A are denoted in what follows 
by < A > n and < A > respectively. 



3 The Chain Approximation Green's Functions 



The GF in the Chain approximation (CHA) are given by [8,12]: 

G&M = r„ ~ w^ ( ":"/ k i^ l2n » (12) 



\Z n £ f,cr) \Zn ^kcr) 


lK(k) 


\ 2 D a 








{z n £/,cr) i,Z n £\ia) 


IK(k) 


?D <7 


- D a V a (k) 






(z n £ f,cr) ( z n &ka) 




?D <7 



GZiZn) = 7 w K \ r^nwr, . (13) 



where 



D a > =< X 00 > + < X a , a , >=R + n fa >. (15) 

and a' = a, a. The lattice slave-boson GF are recovered if we put 
D a = 1 and V k -> rV k = V k in Eqs.(12,13,14). The usual GF for the lattice 
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in the CHA are obtained if we use the bare D a , while D a must be calculated 
self-consistently in the X-boson approach. 



In an earlier paper [15] we considered the atomic limit of the PAM, i.e. with 
a conduction band of zero width, to study both the free energy and N t as a 
function of /x, but here we shall consider the PAM with a wide band. 



4 Calculation of the Lattice Helmholtz Free Energy 

When the total number of electrons Nt, the temperature T and the volume 
V s are kept constant the equilibrium state corresponds to a minimum of the 
Helmholtz free energy, but the same state of equilibrium is obtained by mini- 
mizing the thermodynamic potential Q = — /c B Tln(Q), (where Q is the Grand 
Partition Function) at constant T, V s , and chemical potential ll (this result is 
easily obtained by employing standard thermodynamic techniques). We shall 
then employ Q as the thermodynamic potential that is minimized in the X- 
boson method with Eq. (4) as constraint. The Helmholtz free energy F is then 
given by 



and our first step would be the calculation of f2. 

A convenient way to obtain fl is to employ the method of £ parameter integra- 
tion [15,17,18]. This method introduces a £ dependent Hamiltonian H(£) = 
H + £Hi through a coupling constant £ (with < £ < 1), where H 1 is the 
hybridization in our case. One obtains [12] 



where < A >^ is the ensemble average of an operator A for a system with 
Hamiltonian H(£) and the given values of ll, T, and V s , while Q Q is the ther- 
modynamic potential of the system with £ = 0. This value of £ corresponds to 
a system without hybridization, and one obtains (in the absence of magnetic 
field £ko- = £k and Sf a = if) 



F = N t n + n, 



(16) 




(17) 



fi o = -|E ln t 1 + ex p(-^k)] 



k 



- — ^ In [1 + 2 exp(-/3e»] + N S A(R - 1), 



(18) 



6 



and to calculate £1 in Eq.(17) we use 



(Hi 



2Re 



E V 3*,Ka ( C la X 0a 



Lko- 



(19) 



Employing standard Green's functions techniques [12] we find 

eiK(k)| 2 /J CT 



(Hi)^ — — J duj np(uj) Im 

_™ k,cr 



(lu+ - s fa ) (lu+ - e ka ) - e\VM\ 2 D a 



, (20) 



where uf{x) = 1/ [1 + exp((3x)] is the Fermi-Dirac distribution and u + = 
uj + iO. This equation has an interesting scaling property: it is equal to the 
corresponding expression of the uncorrelated system for the scaled parameter 
K?',k,o- = V Doo-Vjfca (it is enough to remember that by replacing Do a = 1 in 
the GF of the CHA one obtains the corresponding GF of the uncorrelated 
system). Rather than performing the £ and u integrations, we shall use the 
value of the Q u for the uncorrelated system with Vj^ a = \fD^Vj^ a and 
employ Eq. (17) to calculate /„* d£ < H?(g) > 5 = tt u - h'?, where 



Q u = — 

° P 



E In [1 + exp(-/3e k )] + N * In [1 + exp(-/3e / )] 



+ N S A(R - 1) 

(21) 



is the fl u for V^ a = 0. The uncorrelated Hamiltonian for the lattice problem 
is then 



h u = Y. ^A,,c K , + E h flh* + ^ A (^ - 1) 

k,cr j,a 



(22) 



and this Hamiltonian can be easily diagonalized. The corresponding H u can 
be written 



W u = 5><.* a U«i,<r +A(i2-1), 



(23) 



where aj (a ij0 .) are the creation (destruction) operators of the composite par- 
ticles of energies (there are N s + 1 states for each spin a). The calculation 
of Q u is then straightforward [12]: 



n = n + — E In [1 + exp(-/3 cj i)ff )] + A(R - 1), 



(24) 
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where 



TV 

n n = n -Q% = — ^ in 



l + 2exp(-/3e/) 



P "L(l + exp(-/%)) 2 



(25) 



and the eigenvalues of the 7i u are just given by the poles of the GF in the 
CHA (Eq.(12)). In the present case the uj iy(T can be calculated analytically, be- 
cause the Hamiltonian for each spin component er is reduced into N s matrices 
2x2, and one finds 



ha = w±>ff (k) = \ (e Ka + §)) ± l -^e^-e f f + A\V^)\ 2 D a . (26) 



The parameter A is obtained minimizing Q with respect to R [12]. To simplify 
the calculations we shall consider a conduction band with a constant density 
of states, width W = 2D, an hybridization constant V a (k) = V, and e^ c = e^; 
we then obtain 



•'(£ k -E / ) 2 + 4V 2 B„ 



All the correlation effects in Eq. (24) are included in Q D (cf. Eq. (16)), and 
the hybridization redistributes the quasi-particle energies in the same way that 
an uncorrelated system with V = y/D^Vj^ would do; one then expects 
Fermi liquid behavior of the quasi-particles [12]. 



5 Results and Conclusions 



Employing the CHA to calculate the GF of the PAM one discovers an undesir- 
able feature in several regions of the parameter space and for sufficiently low 
T [15]. When the total number of particles in the system N t is plotted against 
the chemical potential \x one finds that there are three possible values of \x for 
each N t in an interval of N t , originating an instability in a small interval of 
/x, characterized by a negative derivative of N t with respect of /i. The problem 
was first observed in a wide conduction band, but further analysis was done 
for a band of zero width [15], because one could then obtain an exact solution 
of the problem and compare with the approximate result. The value of the 
free energy F was then calculated for the three states to analyze their relative 
stability. Here we shall consider the problem for a wide conduction band, and 
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compare the results of the X-boson with those obtained with the CHA and 
with the slave-boson. The main result of the present work is that the X-boson 
does not present the multiple solutions of the CHA. 

In Figure 1 we plot N t (fi) for CHA, slave boson and X-boson for the same 
parameters used in [15]. The multiple solutions of the CHA are apparent, as 
well as the thermodynamically unstable region (part b-c of the Fig. 1, where 
N t decreases when \i increases), i.e. the same behavior that was observed for 
the atomic limit of the PAM [15]. The problem of thermodynamic unstable 
solutions does not appear in the slave boson and X-boson, but the slave boson 
breaks down at higher /j, (// >> ef, magnetic regime) while the X-boson 
presents a continuous evolution in all regimes of the model. 

2.2 i , , , , , 




0.0 



Fig. 1. Total Number of particles N t vs. //, in the three approaches with the following 
parameters: ej = —0.15; W = tt; V = 0.1; T = 0.001. The slave boson breaks down 
when jU f» —0.12. 

In Figure 2 we present the Helmholtz free energy F vs. N t in the slave boson, 
X-boson, and CHA approximations. The CHA present three values of free 
energy F for the same N t that have three values of //, while this problem is 
absent in the X-boson and slave boson methods. 

In Figure 3 we plot the grand thermodynamic potential vs. N t for the slave 
boson and X-boson, at two values of T. In this plot the grand thermodynamic 
potential Q for the slave boson is always lower or equal than for the X-boson, 
and its values increase with temperature in the two methods; its minimum 
values are in the Kondo region at the lowest T. The values of f2 are different 
at the larger value of T for all N t . 

In Figure 4 we show that Q(N t ) has multiple values for the CHA at T = 0.001 
and T = 0.003 when N t ^ 1.5, while that behavior is absent at higher T. Note 
that there is no minimum of Q in the CHA, and that this quantity changes 
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-0.50 1 1 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 

N, 




Fig. 3. Grand thermodynamic potential Q vs. Nt, for the slave boson and X-boson 
approaches with the same parameters of Figure 1, but for two values of the temper- 
ature: T = 0.1 and T = 0.001. 

with N t between maximum and minimum values that are independent of T. 

In Figure 5 we plot N t vs. /x for the CHA at different values of T, and the 
multiple values of /x for a given N t disappear at the larger value of T. 

In Figure 6 we plot p/(/x) (the density of states / on the Fermi level /x) vs. N t 
for the X-boson and slave boson at T = 0.001, the lower value of temperature 
that we consider in Figure 3. From this figure is clear that the Kondo effect 
is present when N t m 1.6 where p/(/x) has a maximum in the two approaches. 
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-0.010 1 1 1 1 1 1 1 

12 3 




Fig. 5. Total Number of particles Nt vs. n, for the CHA approach with the same 
parameters of Figure 1 but at different temperatures values. 

5. 0. 1 Conclusions 

Some years ago we studied the thermodynamic properties of the CHA [15], 
an approximation that was obtained by employing a diagrammatic cumulant 
expansion for the PAM. We then found that for a region of the system param- 
eters at very low temperatures, there are three states with different chemical 
potential /x for the same total number of electrons N t . We observed that one 
of these states is thermodynamically unstable (part b-c of the Fig. 1) be- 
cause N t decreases when /i increases. The X-boson is essentially a generalized 
CHA approximation which satisfies completeness through a minimization of 
the thermodynamic potential as a function of the average occupation of the 
empty state while keeping Eq. (3) as a constraint. As it is shown in Figure (1) 
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1.0 1.2 1.5 1.8 2.0 

N» 

Fig. 6. Density of states / on the Fermi level Pf(n) vs. N t (/j,), for the X-boson 
and Slave boson approaches with the same parameters of Figure 1 . 

this approach does not present thermodynamically unstable states or multiple 
solutions for a given N t . We conclude that the X-boson combines the simplic- 
ity and usefulness of the CHA and slave boson methods, without their more 
obvious defects: the presence of spurious phase transitions of the slave boson 
and the absence of Kondo resonance, of multiple solutions and instabilities, 
as well as the failure to satisfy completeness of the CHA. 
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